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ABSTRACT 

The early evolution of star clusters in the Small Magellanic Cloud (SMC) has been 
the subject of significant recent controversy, particularly regarding the importance and 
length of the earliest, largely mass- independent disruption phase (referred to as "infant 
mortality"). Here, we take a fresh approach to the problem, using an independent, 
homogeneous data set of UBVR imaging observations, from which we obtain the 
SMC's cluster age and mass distributions in a self-consistent manner. We conclude 
that the (optically selected) SMC star cluster population has undergone at most 30 
per cent (Icr) infant mortality between the age range from about (3— 10) Myr, to that 
of approximately (40 — 160) Myr. We rule out a 90 per cent cluster mortality rate per 
decade of age (for the full age range up to 10^ yr) at a > 6(t level. We independently 
affirm this scenario based on the age distribution of the SMC cluster sample. 

Key words: stellar dynamics - globular clusters: general - open clusters and asso- 
ciations: general - Magellanic Clouds - galaxies: star clusters 



1 INTRODUCTION 

The early evolution of the star cluster population in the 
Small Magellanic Cloud (SMC) has been the subject of con- 
siderable recent attention and vigorous debate (e.g., Rafelski 
& Zaritsky 2005; Chandar, Fall & Whitmore 2006; Chiosi et 
al. 2006; Gieles, Lamers & Portegies Zwart 2007). The key 
issue of contention is whether the SMC's star cluster system 
has been subject to the significant early cluster disruption 
processes observed in "normal", interacting and starburst 
galaxies commonly referred to as "infant mortality" and "in- 
fant weight loss" . Chandar et al. (2006) argue that the SMC 
has been losing up to 90 per cent of its star clusters per 
decade of age, at least for ages from ~ 10^ up to ~ 10^ 
yr, whereas Gieles et al. (2007) conclude that there is no 
such evidence for a rapid decline in the cluster population, 
and that the decreasing number of clusters with increas- 
ing age is simply caused by fading of their stellar popula- 
tions. They contend that the difference between their results 
was due to Chandar et al. (2006) assuming that they were 
dealing with a mass-limited sample, whereas it is actually 
magnitude- limited. In fact, this is not entirely correct; Chan- 
dar et al. (2006) analyse the full magnitude-limited sample 
and conclude that it is approximately surface-brightness lim- 
ited. They then compare the cluster age distribution of the 
full sample (expressed in units of AN^i/At, i.e., the number 
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of clusters per unit time period) to that of a subsample for 
masses > 10'^ M0 (which they do not analyse in the same 
manner), and suggest both to be similar, although the latter 
is much flatter^, hence giving rise to the discrepancy between 
their results and those of Gieles et al. (2007). Both stud- 
ies are based on the same data set, the Magellanic Clouds 
Photometric Survey (MCPS; Zaritsky, Harris & Thompson 
1997). 

The main contribution of this paper to the ongoing de- 
bate is two-fold: (i) We revisit the SMC's early star clus- 
ter evolution using an alternative approach; and (ii) we use 
independently obtained ages and masses based on an inde- 
pendent data set, i.e., the UBVR photometric survey of the 
Magellanic Clouds by Massey (2002), originally analysed by 
Hunter et al. (2003). We conclude that there is indeed only 
marginal evidence for infant mortality in the SMC star clus- 
ter sample, supporting the careful analysis of Gieles et al. 
(2007). In Sect. 2 we first briefly introduce the concept of 
cluster infant mortality. We discuss our observational data 
and the basic analysis leading to the age and mass estimates 
in Sect. 3. In Sect. 4 we justify our choice of age ranges to 
construct cluster mass functions (CMFs). Finally, in Sect. 



^ Although Chandar et al. (2006) suggest that their sample is 
roughly mass limited, they also note that the mass-limited sub- 
sample, constrained to clusters with masses log(M(.i/M0) > 3.5, 
shows a flatter age distribution. 
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5, we present our case for the absence of significant cluster 
infant mortality. 



2 CLUSTER INFANT MORTALITY 

Observations of increasing numbers of interacting and star- 
burst galaxies, including the Antennae system, M51 and 
NGC 3310, show a significantly larger number of young 
(< 10 — 30 Myr) star clusters than oxpoctod from a simple 
extrapolation of the cluster numbers at older ages, taking 
into account the observational completeness limits and the 
effects of sample binning, and under the additional, simplify- 
ing assumption that the star cluster formation rate (CFR) 
has been roughly constant over the host galaxy's history 
(e.g., dc Grijs et al. 2003b; Whitmore 2004; Bastian et al. 
2005; Fall, Chandar & Whitmore 2005; Mengel et al. 2005; 
Chandar et al. 2006; see also Whitmore, Chandar & Fall 
2007 for a presentation of earlier results, and de Grijs & 
Parmcnticr 2007 for a review). This significant ovordcnsity 
remains, even in view of the presence of a recent burst of 
star cluster formation in many of these galaxies. 

These observations have prompted a flurry of activity 
in the area of star cluster disruption processes. This has led 
to suggestions that cluster systems appear to be affected by 
a disruption mechanism that acts on very short time-scales 
(< 10 — 30 Myr) and which may be mass-independent - 
at least for masses in excess of ~ 10* M© (e.g.. Fall et al. 
2005; Bastian et al. 2005; Fall 2006). This fast disruption 
mechanism, which is thought to effectively remove axound 
50 (Goodwin & Bastian 2006; although their sample is very 
probably biased), 70 (Bastian et al. 2005; Mengel et al. 2005) 
or even 90 per cent (Lada & Lada 1991; Whitmore 2004; 
Whitmore et al. 2007) of the youngest clusters from a given 
cluster population, is thought to be the rapid removal of the 
intracluster gas on a time-scale of ~ 5 Myr, the signatures of 
which have been seen in several clusters (Bastian & Goodwin 
2006). The observational effect resulting from this rapid gas 
removal has been coined cluster "infant mortality" (Lada & 
Lada 2003); it was originally reported in the context of the 
number of very young embedded clusters, compared to their 
older, largely gas-free counterparts in the Milky Way. 

The general consensus emerging from recent studies is 
that rapid gas removal from young star clusters is likely to 
leave the clusters super-virial and hence lead to the rapid 
disruption of many clusters (see, e.g., Goodwin 1997a,b; 
Boily & Kroupa 2003a,b; Goodwin & Bastian 2006; see also 
de Grijs & Parmentier 2007 for a review). This leaves sur- 
viving clusters more susceptible to destruction (Vesperini & 
Zepf 2003; Bastian et al. 2005; Fall et al. 2005). 

As described by Goodwin & Bastian (2006, and refer- 
ences therein) the effect of gas removal is to rapidly decrease 
the potential well in which the stars reside. The cluster will 
expand in an attempt to return to virial equilibrium. If the 
virial ratio of the stars after gas expulsion is > 3 the clus- 
ter will be unable to return to an cquilibirum and will be 
destroyed^. Clusters with a higher effective star- formation 

This is usually given in terms of a star-formation efficiency of 
~ 30 per cent. However, as noted by Goodwin & Bastian (2006), it 
is the virial ratio of the stars after gas expulsion that is the crucial 
parameter, and this can only be related to the star-formation 



efficiency than around 30 per cent will survive, but may un- 
dergo significant "infant weightloss" (loosing in excess of 50 
per cent of their initial mass in some cases). The signature 
of infant weightloss has been observed in several young clus- 
ters (Bastian & Goodwin 2006). The time-scale over which 
a cluster will be destroyed, or attain a new (lower-mass) 
equilibrium configuration is 10 — 40 Myr (depending on the 
effective star- formation efficiency and the cluster mass). 

The 10 — 40 Myr time-scale of gas removal-induced in- 
fant mortality and infant weightloss has important conse- 
quences for the analysis and interpretation of the data in 
this paper. Clusters undergoing expansion will have decreas- 
ing surface brightnesses, thus reducing their chances of be- 
ing detected as they grow older. However, some clusters will 
recollapse after 10 — 40 Myr, which may bring them back 
into the sample. In addition, the speed at which clusters 
are lost from the sample would be expected to depend on 
their (initial) mass. Lower-mass clusters which arc initially 
only just above the detection limit will drop out of the sam- 
ple very quickly, whilst larger clusters may remain in the 
sample (albeit with a lower surface brightness) for longer. 
It is also almost impossible (without extensive observations 
of the surface brightness profiles and/or dynamical state of 
the cluster) to determine which clusters that are present 
in this age range will survive gas expulsion, and which are 
headed for destruction. Thus, the interpretation of the num- 
bers and mass function of clusters in the age range 10 — 40 
Myr is fraught with problems (in addition to these "phys- 
ical" problems, age determinations for clusters in this age 
range also cause problems; see below). 



3 A HOMOGENEOUS PHOTOMETRIC 
DATABASE 

The basis for our detailed re-analysis of the SMC star clus- 
ter system is provided by the UBVR broad-band spectral 
energy distributions (SEDs) of Hunter et al. (2003), based 
on Masscy's (2002) CCD survey of the MageUanic Clouds. 

In a scries of recent papers, we developed a sophis- 
ticated tool for star cluster analysis based on broad-band 
SEDs, Analysed, which we tested extensively both inter- 
nally (dc Grijs et al. 2003a,b; Anders et al. 2004) and ex- 
ternally (de Grijs et al. 2005), using both theoretical and 
observed young to intermediate-age (< 3 x 10® yr) star clus- 
ter SEDs, and the galev "simple" stellar population (SSP) 
models (Kurth et al. 1999; Schulz et al. 2002). The accu- 
racy has been further increased for younger ages by the 
inclusion of an extensive set of nebular emission lines, as 
well as gaseous continuum emission (Anders & Fritze-v. Al- 
vensleben 2003). Wc concluded that the relative ages and 
masses within a given cluster system can be determined to 
a very high accuracy, depending on the specific combination 
of passbands used (Anders et al. 2004). Even when com- 
paring the results of different groups using the same data 
set, wc can retrieve any prominent features in the cluster 
age and mass distributions to within A(log(Age/yr)) < 0.35 
and A(log(Mci/M0)) < 0.14, respectively (de Grijs et al. 

efficiency in a simple way if the stars and gas were initially in 
virial equilibrium with one another. For this reason, Goodwin & 
Bastian (2006) use the term "effective star formation efficiency" . 
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Figure 1. Distribution of the SMC clusters in the log(age) vs. 
log(mass) plane. Overplotted is the expected detection limit based 
on stellar population synthesis for a 50 per cent completeness 
limit of My = —4.5 mag, assuming no extinction. For a nomi- 
nal extinction of Ay = 0.08 mag (assuming the Calzetti atten- 
uation law), the detection limit will shift to higher masses by 
A log(Mt.i/MQ) = 0.03, which is well within the uncertainties as- 
sociated with our mass determinations (see de Grijs & Anders 
2006). The features around 10 Myr are caused by the appearance 
of red supergiants in the models. The age limits used to gener- 
ate the different panels in Fig. 2 are shown as the vertical dash- 
dotted lines; the various subsets are also cross-linked between the 
figures using the panel indications from Fig. 2. The horizontal 
dash-dotted lines indicate the 50 per cent completeness limits in 
mass for each of the age-selected subsamples. 



2005), which confirms that we understand the uncertainties 
associated with the use of our ANALYSED tool to a very 
high degree. 

In de Grijs & Anders (2006) we presented newly and ho- 
mogeneously redetermined age and mass estimates for the 
entire Large Magellanic Cloud (LMC) star cluster sample 
covered by the Massey (2002) data. Based on the comparison 
of our results in de Grijs & Anders (2006) with those pub- 
lished previously in a range of independent studies (mostly 
based on spectroscopic or isochrone analyses) , and addition- 
ally on a detailed assessment of the age-metallicity and age- 
extinction degeneracies, we concluded that our broad-band 
SED fits yield reliable ages, with statistical absolute uncer- 
tainties within Alog(Age/yr) ~ 0.4 overall. Here, we extend 
this to the SMC cluster sample, using the same age-dating 
technique as described above. 

Our cluster age and mass determinations assume an av- 
erage metallicity of Z — 0.008 (where Zq = 0.020), and a 
mean foreground extinction E{B — V) = 0.08 mag. We will 
justify both of these choices below. In Fig. 1 we show the 
distribution of our SMC cluster sample in the log(age) vs. 
log(mass) plane; the adopted 50 per cent completeness limit 
is overplotted. We have also indicated the regions in this 



plane from which we have drawn statistically complete sub- 
samples, which we will discuss in detail in Sect. 4. 

The determination of the 50 per cent completeness limit 
of the SMC cluster data is in essence based on a close inspec- 
tion of the cluster photometry contained in Hunter et al.'s 
(2003) fig. 11. These authors selected their sample from the 
catalogue of Pietrzynski et al. (1998), matched to the ob- 
servational field of view of the Massey (2002) data. There- 
fore, our completeness is that of this catalogue; Hunter and 
her team did not quantify the completeness levels them- 
selves (D. Hunter, priv. comm.), although they discuss an 
observed fading limit. However, for our analysis it is impor- 
tant to understand the sample incompleteness affecting our 
observations. As such, we adopted the conservative approach 
that the present-day SMC cluster luminosity function (CLE; 
see Hunter et al.'s fig. 11) is best represented by a power- 
law function in luminosity. Based on this assumption, we 
used the same observational data as used by Hunter et al. 
(2003) to determine the 50 per cent completeness limit at 
AIv ~ —4.5 ± 0.2 mag (based on a power-law fit to the clus- 
ters brighter than My ~ —5 mag; varying this lower limit 
by a few tenths of a magnitude does not result in a signifi- 
cant change), i.e., at the same level as Hunter et al.'s (2003) 
observed fading limit. We note that if the underlying CLE 
is not a single power law over the entire observed luminosity 
range, the limit we adopt following this approach is in fact 
a lower limit. In the latter case the observations will likely 
be more complete than estimated here. Since the adopted 
50 per cent completeness limit describes the lower envelope 
of the distribution of our SMC cluster sample very well, we 
are confident that our approach is reasonable. In addition, 
we point out that a variation in the magnitude limit of 0.2- 
0.4 mag will shift the mass limit by at most 0.1-0.2 dex, 
which clearly is still within our range of uncertainties. The 
magnitude of the shift expected when going from the 50 to 
the 90 per cent completeness limit in the SMC disc is of the 
same order, ~ 0.5 mag. Finally, we point out that it is most 
likely that the completeness of our SMC cluster sample is in 
fact determined by the (7-band observations. Erom a direct 
comparison of the [/-band and the \/-band data, we derive a 
50 per cent completeness in the U band at Mu = —5.0 ±0.3 
mag. 

Chiosi et al. (2006) recently analysed the star-formation 
history in the SMC in detail using an independently selected 
star cluster sample. Where possible, they derive the extinc- 
tion towards individual clusters based on colour-magnitude 
diagram analysis, and otherwise assume a mean extinction 
E{B-V) = 0.08 mag, following Tumlinson et al. (2002) and 
Hunter et al. (2003; see also Rafelski & Zaritsky 2005). We 
have adopted the same average extinction value to our SMC 
cluster sample, using the Calzetti attenuation law (Calzetti 
1997, 2001; Calzetti et al. 2000; Leitherer et al. 2002) with 
Ry = 4.05. 

Rafelski & Zaritsky (2005) obtained SMC cluster ages 
of a small cluster sample on the basis of three sets of mod- 
els, for metallicities oi Z = 0.001,0.004 and 0.008. They 
concluded that some of the lowest-metallicity models could 
be rejected and adopted Z — 0.008 as an appropriate mean 
metallicity for their SMC cluster sample. Chiosi et al. (2006) 
also adopted this metallicity, but for their younger sample 
clusters. For the older (> 1 — 2 Gyr) clusters, they assumed 
Z = 0.004, as was also done by Hunter et al. (2003). How- 
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ever, as shown in Fig. 1, the large majority of our sample 
clusters (and in particular the subpopulations we will anal- 
yse in more detail below) are younger than ~ 1 Gyr, so that 
we adopt Z = 0.008 as the mean metallicity for our SMC 
cluster sample. 



4 THE CLUSTER MASS FUNCTION 

In dc Grijs & Anders (2006) we presented the cumula- 
tive CMFs of the LMC star clusters younger than cer- 
tain age limits. We concluded that while the older cluster 
(sub) samples arc characterised by CMF slopes consistent 
with the Q ~ 2 slopes generally observed in young star clus- 
ter systems - where a is defined as iV(Mci) oc M^^" ~ the 
youngest mass and luminosity-limited LMC cluster subsets 
show shallower slopes (at least below masses of a few x 10^ 
M0). We noted that we could not disentangle the unbound 
from the bound clusters at the youngest ages. This is what 
we set out to do here for the SMC cluster system. 

In order to achieve this goal, we present the CMFs for 
subsets of our SMC cluster sample in Fig. 2, where the clus- 
ter subsamples were selected based on their age distribu- 
tions. The age limits used to generate the different panels in 
Fig. 2 are shown as the vertical dash-dotted lines in Fig. 1; 
the various subsets are also cross-linked between the figures. 

A closer look at Fig. 1 reveals that, because of the vari- 
ation in the observational detection limit as a function of 
age, the lower-mass limits of our cluster subsamples differ. 
Thus, the CMFs presented in Fig. 2 are statistically com- 
plete above different mass limits, as indicated by the hor- 
izontal dash-dotted linos in Fig. 1 and the vertical dotted 
lines in Fig. 2. 

In all panels of Fig. 2, we have overplotted CMFs with 
the canonical slope of a = 2 (corresponding to a slope of —1 
in units of d log(Mci/M0)/d log(A'ci), used in those panels). 
We have only shifted and scaled these lines vertically, as 
justified below. 

We emphasise that for the star cluster infant mortal- 
ity study performed here, we need to choose the age ranges 
of our cluster subsamples carefully, for both physical rea- 
sons and also because of the discrete nature of the model 
isochrones. Regarding the latter, it is well known that broad- 
band SED fitting results in artefacts in the cluster age distri- 
bution. This is predominantly caused by specific features in 
the SSP models, such as the onset and presence of red giant 
branch or asymptotic giant branch (AGB) stars at, respec- 
tively, ~ 10 and ~ 100 Myr (e.g., Bastian et al. 2005). Alter- 
natively, both the agc-metallicity and the ago-oxtinction de- 
generacies will affect the resulting cluster age distributions, 
thus also leading to artefacts in the data (e.g., de Grijs et 
al. 2003b; Anders et al. 2004). We have attempted to avoid 
placing our age range boundaries around ages (and, where 
possible, have taken account of the uncertainties in ago in 
doing so) where the effects of such artefacts might seriously 
impede the interpretation of the results. For instance, one 
can see a clear artefact in the cluster age distribution (which 
we will refer to as a "chimney") at log(Age/yr) ~ 7.2(~ 16 
Myr); the average uncertainties for those ages arc of order 
a few Myr, so that we decided to limit our youngest cluster 
subsample to clusters younger than 10 Myr. If, instead, we 
had adopted an age limit at log(Age/yr) = 7.17(15 Myr), we 



would have had marginally better statistics, but our analysis 
would be affected by the unknown effects of the age uncer- 
tainties associated with this chimney (see Goodwin et al., in 
prep., for a detailed discussion of the issues involved). 

The rationale for adopting as our youngest subsample 
all clusters with ages < 10 Myr is that at these young ages, 
the vast majority of the star clusters present will still be 
detectable, even in the presence of early gas expulsion (e.g., 
Goodwin & Bastian 2006) - as long as they are optically 
conspicuous. The CMF of this subsample is shown in Fig. 
2a. 

Fig. 2c includes our sample clusters with ages in excess 

of 40 Myr, up to 160 Myr. While the upper age limit ensures 
the full inclusion of the clusters affected by the onset of the 
AGB stage, its exact value is rather unimportant for our 
analysis, and it was mainly determined by the need to have 
reasonable statistics in this (and the upper) age range in Fig. 
2d. The lower age limit of this subsample is crucial, however. 
As shown by Goodwin & Bastian (2006), most dissolving 
clusters will have dispersed by an age of ^ 30 Myr, while the 
surviving clusters will have returned to an equilibrium state 
by ~ 40 Myr, when some of the early expansion will have 
been reversed, depending on the effective star-formation ef- 
ficiency. This latter age is therefore a good lower boundary 
to assess the surviving star cluster population. 

We explicitly exclude any star clusters aged between 10 
and 40 Myr from our analysis. In this age range, which is 
shown in Fig. 2b for completeness, it is likely that dissolv- 
ing star clusters that will not survive beyond about 30-40 
Myr might still be detectable and therefore possibly contam- 
inate our sample. In addition, this is the age range in which 
early gas expulsion causes rapid cluster expansion, before 
settling back into equilibrium at smaller radii; because of 
the expanded nature of at least part of the cluster sample, 
we might not be able to detect some of the lower-luminosity 
(and hence lower-mass) clusters that would again show up 
beyond an age of ^ 40 Myr. At the same time, the effects of 
"infant weightloss" (Weidner et al. 2007) will further confuse 
the analysis in this age range (see Section 2 for details). 

5 IS THERE EVIDENCE FOR CLUSTER 
INFANT MORTALITY IN THE SMC? 

5.1 Young and intermediate-age clusters 

In Fig. 2a, we have included the best-fitting CMF slope 
(dash-dotted line) , in addition to the canonical a = 2 CMF 
slope (dashed line). Both slopes are the same, within the 
uncertainties. This also shows that the SMC's CMF at the 
youngest ages is consistent with an a = 2 power law down 
to cluster masses of ~ 125 M0, within the (Poissonian) un- 
certainties. 

In the simplest case, in which the cluster formation rate 
has remained roughly constant throughout the SMC's evo- 
lution (see, e.g., Boutloukos & Lamers 2003, their fig. 10; see 
also Gieles et al. 2007), the number of clusters would simply 
scale with the age range covered. In Figs. 2b, c and d, we 
show the canonical a = 2 CMF scaled from the best-fitting 
locus iu Fig. 2a by the difference in (linear) age range be- 
tween the panels. The main uncertainties introduced by this 
method are (i) fluctuations caused by small-number statis- 
tics in the youngest age range (the effects of which will be 
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Figure 2. CMFs for statistically complete SMC cluster subsamples. Age and mass ranges are indicated in the panel legends; the vertical 
dotted lines indicate the lower mass (50 per cent completeness) limits adopted. Error bars represent simple Poissonian errors, while the 
dashed lines represent CMFs of slope a = 2, shifted vertically as described in the text. The dash-dotted lines represent the best-fit CMFs 
over the relevant mass range. The panel indicators refer to Fig. 1. 



enhanced when scaling the young-age CMF to a greater age 
range), and (ii) the exact length of the youngest age range 
(especially considering the necessarily short extent of our 
youngest age bin). While our GALEV SSP models start at an 
age of 4 Myr, the actual ages of a small subset of our sample 
clusters might be somewhat younger. This introduces an ar- 
tificial concentration of clusters at our youngest model age, 
as can be seen in Fig. 1. It is, unfortunately, not straightfor- 
ward to remedy this situation based on broad-band imaging 
observations alone. However, we note that it may take up 
to 3-5 Myr for an embedded cluster to clear a cavity in 
its natal gas cloud for its stars to become visible at optical 
wavelengths (see, e.g., Plante & Sauvage 2002 for a review of 
embedded young massive star cluster observations). There- 
fore, we adopt the conservative working assumption that our 
youngest age range runs from 3-10 Myr. The scaling from 
the youngest age bin to that covering [40,160] Myr (Fig. 2c, 
i.e., the most important age range for our CMF comparison 
in the context of our infant mortality analysis) is therefore 
a factor of ~ 17, or Alog(iVci) = 1.234. 

Despite the caveat regarding the absence of embedded 
star clusters in our youngest subsample, we argue that this 
has a negligible effect on the CMF presented in Fig. 2a, 
because of their very small number. Recent, homogeneous 



observations of the SMC using the Spitzer Space Telescope 
in a number of mid-infrared passbands (Bolatto et al. 2007) 
have shown that the vast majority of the embedded sources 
are low-mass (<C 100 Mq) young stellar objects rather than 
more massive clusters and associations (we point out that 
for the comparison done here, we are mostly interested in 
clusters with masses greater than 10^ M©). The possible 
exceptions to this rule are few, and include the four youngest 
~ 3 Myr-old SMC clusters, NGC 299, NGC 346, NGC 376, 
and NGC 602 (e.g., Sabbi et al. 2007). 

The scaled canonical CMF in Fig. 2c is an almost per- 
fect fit to the observational CMF. Although the best-fitting 
CMF slope is dlog(Mc:i/MQ)/dlog(iVc:i) = -0.82±0.18, this 
compares to dlog(Mci/M0)/dlog(iVci) = -1.01±0.20 if we 
ignore the lowest-mass clusters at log(Mci/M0) < 3.2, where 
there may be residual incompleteness effects (see the selec- 
tion area in Fig. 1 compared to the age-dependent detection 
limit). 

This very good match between the observed CMF for 
the age range from 40-160 Myr (Fig. 2c) and the scaled 
CMF from Fig. 2a implies that the SMC cluster system has 
not been affected by any significant amount of cluster infant 
mortality for cluster masses greater than a few x 10^ Mq . 
Based on a detailed assessment of the uncertainties in both 
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the CMFs and the age range covered by our youngest sub- 
sample, we can Umit the extent of infant mortality between 
the youngest and the intermediate age range to a maximum 
of < 30 per cent (1(t). We rule out a 90 per cent (infant) 
mortality rate per decade of age at a > 6(J level. This result 
is in excellent agreement with that of Gieles et al. (2007) ; it 
is, however, in direct contradiction to the claim of Chandar 
et al. (2006) that the SMC cluster system has undergone 
sustained destruction at very high rates (up to 90 per cent 
per decade in logarithmic age) for the full age range up to 
~ 1 Gyr, although we note that Chandar et al. (2006) do 
not include the youngest SMC clusters in their analysis. 

As an important caveat, we remind the reader that the 
main underlying assumption leading to this result is the no- 
tion that the SMC's CFR has been approximately constant 
over the time-scale of ~ 1 Gyr. If this were seriously in error, 
in order for this to give rise to the result reported here, the 
SMC's average CFR must have been significantly enhanced 
in the 40 — 160 Myr-old age range, by up to an order of 
magnitude, compared to that at present. There is no clear 
evidence, in either the current data set or the earlier work by 
Boutloukos & Lamers (2003; see also Gieles et al. 2007), nor 
in the age distribution of the field stars (Chandar et al. 2006; 
Chiosi & Vallenari 2007) , to suggest that this is the case. In 
fact, if anything, we might expect an enhanced CFR around 
the time of the last close encounter between the SMC and 
the LMC, some 200-500 Myr ago (see, e.g.. Heller & Rohlfs 
1994; Gardiner & Noguchi 1996; see also Chiosi et al. 2006, 
but sec Chiosi & Vallenari 2007 for an alternative interpreta- 
tion), i.e., significantly longer ago than the age range probed 
by our intermediate-age clusters in Fig. 2b. 

For completeness, we also show the best-fit power-law 
CMF, as well as the scaled canonical CMF, in Fig. 2b for 
the age range between 10 and 40 Myr. Although we strongly 
caution against placing too much significance on the analy- 
sis of the CMF in this age range, for the reasons outlined in 
Section 4, it is interesting to note that the scaled canonical 
CMF does in fact seem to describe the extremities of this 
CMF reasonably well. However, at intermediate masses (a 
few xlO''* - lO"* Mq) the observed CMF exceeds the theo- 
retical prediction for a constant cluster formation rate. Al- 
though the effects of cluster expansion and infant weightloss 
most likely contribute to confusing the emerging picture, the 
main cause of this discrepancy is owing to the artificial chim- 
ney at log(Age/yr) ~ 7.2. The cluster numbers in this age 
range are dominated by this artefact, so it is important that 
we understand in which sense this affects our results. Be- 
cause of the discreteness of the isochrones in our SSP models 
around this age, and the tendency for the broad-band fitting 
routine to iterate to a local minimum solution, most (but 
not all) of the clusters in this chimney should have been as- 
signed somewhat greater ages. Because their ages have been 
underestimated (by up to 0.1 0.2 dex in logarithmic age), 
the associated masses have also been underestimated, by 
a similar amount. Unfortunately, until more detailed SSP 
models become available, there is no easy way out. How- 
ever, a qualitative assessment suggests that if we were able 
to correct for this chimney, the derived masses of at least a 
fraction of the clusters affected would be greater, and thus 
that the apparent excess in Fig. 2b would be redistributed 
towards greater masses. The result would be a smoother 



CMF, more akin to the scaled canonical a = 2 CMF of Fig. 
2a. 

The alternative interpretation, i.e., that the star clus- 
ter formation rate in the SMC has undergone a significant 
increase in the age range between 10 and 40 Myr appears 
to be effectively ruled out by the complementary analysis 
of Chiosi et al. (2006). In fact, these authors find evidence 
suggesting the contrary, i.e., that the SMC cluster popular 
tion has seen periods of enhancement during the first ~ 15 
Myr, and at around 90 Myr - this implies that in the age 
range of interest here, the cluster formation rate found by 
Chiosi et al. (2006) was in fact reduced with respect to the 
most recent cluster forming episode (roughly equivalent to 
the cluster sample shown in Fig. 2a). Equivalcntly, neither 
Rafelski & Zaritsky (2005), nor either Chandar et al. (2006) 
or Gieles et al. (2007) find an episode of increased cluster 
formation at a few x lO'^ yr, despite the fundamentally dif- 
ferent conclusions drawn by the authors of the latter two 
studies. 



5.2 The oldest sample clusters 

Finally, in Fig. 2d we show the combined SMC CMF 

for clusters from 160 Myr up to 1.0 Gyr, as well as the 
scaled canonical CMF. The latter matches the highest- 
mass (log(Mci/M0) >4.8) part of the observed CMF, 
but significantly overpredicts the number of lower-mass 
(log(Mci/M0) < 4.8) clusters. This fiattening of the CMF 
with respect to our intermediate age range (Fig. 2c) ev- 
idences the increased importance of mass-dependent clus- 
ter disruption, as described in detail by, e.g., Lamers et al. 
(2005; see specifically their fig. 4). However, we note that 
the calculations of Lamers et al. (2005) are based on co- 
eval stellar populations, while our oldest age bin contains a 
mixture of differently aged clusters. Nevertheless, it is ap- 
parent from fig. 4 of Lamers et al. (2005) that the flattening 
of the CMF increases as a function of age, starting at early 
ages. In a mixed intermediate-age population, the individ- 
ual roughly coeval subpopulations giving rise to the overall 
CMF therefore lead to a flattened CMF with respect to the 
initial cluster mass distribution (which we have shown in 
the SMC to be equivalent to the canonical a = 2 power- 
law CMF; see Fig. 2a). A comparison of our observations 
with the Lamers et al. (2005) calculations is therefore rel- 
evant to flrst order. A similar flattening of the CMF with 
increasing age is also predicted as owing to the effects of 
the underlying galactic tidal field (e.g., Vesperini & Heggie 
1997; Baumgardt & Makino 2003). We will discuss these 
older clusters in more detail in a follow-up paper (Goodwin 
et al., in prep.), in which we will discuss the entire evolution- 
ary sequence of the star cluster systems of both Magellanic 
Clouds, and where we aim to understand the ph/ysics driving 
the early evolution of these star clusters systems. 

5.3 Additional supporting evidence for little early- 
disruption 

In Fig. 3 we show the SMC cluster age distribution expressed 
in units of dNd/dt, i.e., the number of cluster per unit time- 
scale (for which we adopt 10® yr). To first order, our age dis- 
tribution is similar to that based on the Rafelski & Zaritsky 
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Figure 3. Age distribution of the SMC cluster sample in units of 
cluster number per Myr. We show four different (sub)samples, in- 
cluding the full magnitude-limited SMC sample, and three mass- 
limited subsamples, as indicated in the figure legend. The mass- 
limited subsamples have been shifted vertically by the constant 
offsets indicated on the right-hand side for reasons of clarity 
(without these offsets, the data points would all overlap at the 
oldest ages). For that same reason, we have also connected the 
data points for each of the (sub)samples. The mass-limited sub- 
samples are > 50 per cent complete to the left of the vertical 
dashed lines at the bottom of the figure, where the numbers refer 
to the 50 per cent completeness limits for a given range, expressed 
in log{M(.i/MQ). The vertical error bars are simple Poissonian er- 
rors; the horizontal error bars indicated the age range used for the 
generation of these data points. Finally, the dashed arrow shows 
the expected effects due to fading of a cluster sample made up 
of SSPs, based on the GALEV SSP models (see also Gieles et al. 
2007) ; the dash-dotted arrow represents the combined effects of a 
fading cluster population and 90 per cent cluster disruption per 
decade in log(Age yr~^), up to ages of 1 Gyr. 



(2005) data used by both Chandar et al. (2006) and Gieles 
et al. (2007). It is also roughly similar to the age distribution 
derived independently by Chiosi et al. (2006). 

Gieles et al. (2007) analysed their cluster age distri- 
bution very carefully, and found little evidence for infant 
mortality in the SMC cluster system. They showed that the 
decline in dNc\/dt in their sample could be attributed en- 
tirely to evolutionary fading of the cluster population, irre- 
spective of which SSP models are used. They derived that 
for a magnitude- limited sample, as also discussed in this pa- 
per, the decline in dNd/dt as a function of age is graphically 
described by a slope of —0.72, if the cluster ages are based 
on the GALEV SSP models. In Fig. 3, we show the expected 
effects of evolutionary fading of the cluster population as the 
dashed arrow. It is immediately clear that the decline in the 
age distribution up to log(Age yr~^) ~ 7.8 can indeed be 
entirely attributed to evolutionary fading. This is supported 
by the mass- limited subsamples shown in Fig. 3: for all mass 
ranges, they show an essentially flat age distribution up to 



~ 10 yr. This is in support of the results of both Gieles et 
al. (2007), and Chandar et al. (2006; their fig. 1), although 
the latter authors favoured a different interpretation. 

The age distribution at log(Age yr~^) ~ 8.2 ± 0.3 (cov- 
ering the age range from about 80 to 320 Myr) falls below 
the fading line, however. For a constant cluster formation 
rate over this entire period, and if we normalise the age dis- 
tribution at our youngest data point, we would need the 
SMC cluster sample to have suffered from ~ 20 — 50 per 
cent disruption in order to match the observations. We can 
firmly rule out a constant ~ 90 per cent disruption rate per 
decade in age, up to an age of 1 Gyr. The expected effects of 
evolutionary fading combined with a 90 per cent disruption 
rate are shown as the dash-dotted arrow in Fig. 3. The ar- 
row clearly does not fit the observed age distribution, if we 
require it to pass through our youngest data point. We note, 
however, that the slope of this latter arrow is very similar 
to that of the age distribution of the full SMC sample for 
ages in excess of a few x 10* yr, when secular disruption is 
likely to take over. We will discuss these older age ranges in 
detail in Goodwin et al. (in prep.). 

In summary, we set out to shed light on the contro- 
versy surrounding the early evolution and disruption of star 
clusters in the SMC. We embarked on a fresh approach to 
the problem, using an independent, homogeneous data set 
of UBVR imaging observations, from which we obtained 
the cluster age distribution in a self-consistent manner. We 
conclude that the optically selected SMC star cluster pop- 
ulation has undergone at most ~ 30 per cent (Icr) "infant 
mortality". Using the age distribution of the SMC cluster 
sample in units of the number of clusters observed per unit 
time-scale, we independently confirm this scenario. Gieles et 
al. (2007) reached a similar conclusion. 
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